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ABSTRACT 

We report our discovery and early time optical, near-infrared, and radio 
wavelength follow-up observations of the afterglow of the gamma-ray burst 
GRB 021211. Our optical observations, beginning 21 min after the burst trigger, 
demonstrate that the early afterglow of this burst is roughly three magnitudes 
fainter than the afterglow of GRB 990123 at similar epochs, and fainter than al- 
most all known afterglows at an epoch of 1 d after the GRB. Our near-infrared 
and optical observations indicate that this is not due to extinction. Combining 
our observations with data reported by other groups, we identify the signature 
of a reverse shock. This reverse shock is not detected to a 3-0" limit of 110 /zJy 
in an 8.46-GHz VLA observation at t = 0.10 d, implying either that the Lorentz 
factor of the burst 7 < 200, or that synchrotron self-absorption effects dominate 
the radio emission at this time. Our early optical observations, near the peak of 
the optical afterglow (forward shock), allow us to characterize the afterglow in 
detail. Comparing our model to flux upper limits from the VLA at later times, 
t > 1 week, we find that the late-time radio flux is suppressed by a factor of two 
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relative to the >80//Jy peak flux at optical wavelengths. This suppression is 
not likely to be due to synchrotron self-absorption or an early jet break, and we 
suggest instead that the burst may have suffered substantial radiative corrections. 

Subject headings: galaxies: high-redshift - gamma rays: bursts 

1. Introduction 

Gamma-ray burst (GRB) afterglow emission on timescales of hours to days, the his- 
torical standard for ground-based follow-up observations, can provide a robust measure of 
the global parameters of the GRB such as its total explosive yield and the density of the 
circumburst medium on length scales of r > 10^^ cm (Sari et al. 1998; Wijers & Galama 
1999). 

Early emission, on a timescale of tens of minutes and less, can offer insight into the 
details of the explosion, including the relativistic Lorentz factor of the burst ejecta (Sari & 
Piran 1999b) and the distribution of the circumburst medium on length scales r > 10^^ cm. 
If the GRB progenitor is a massive star, then this regime is expected to show the clear 
signature of a dense stellar wind (density p oc r~^; Chevaher & Li 2000). 

These two unique diagnostics, among others, have motivated searches for early optical 
emission from GRBs. The first success came with detection of early optical (Akcrlof ct al. 
1999) and radio (Kulkarni et al. 1999) emission from GRB 990123 (Briggs et al. 1999). This 
emission has been interpreted in terms of a reverse shock plowing back through the GRB 
ejecta (Sari & Piran 1999a; Meszaros & Rees 1999). 

After nearly four years, early emission was recently discovered again from GRB 021004 
(Fox et al. 2003) and GRB 021211 (Fox & Price 2002). These discoveries were made possible 
by the prompt alerts of the HETE-2 saiellite (Ricker et al. 2002), and the wide field and high 
sensitivity of the 48-inch Palomar Oschin Schmidt telescope (P48) equipped with the NEAT 
camera (Pravdo et al. 1999). The increase in the number of robotic telescopes has allowed for 
intensive follow-up of both of these bursts during the critical early minutes to hours. Here, 
we report the discovery of early optical emission and subsequent multi-wavelength follow-up 
observations of GRB 021211. 
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2. Observations 

GRB 021211 (HETE Trigger #2493) was detected at UT 11:18:34 on 11 Dec 2002. 
With a duration of 6 s and a fluence of lO^^erg cm^^ (8-40 keV), the event is a typical 
long duration GRB (Crew et al. 2002). At the time of electronic notification the HETE-2 
localization, although large (30-arcminute radius), was covered by the 1.1° x 1.1° field- 
of-view of NEAT. We identified a new source that was not present in the Palomar Sky 
Survey (Fig. 1) and announced this as a possible optical counterpart of GRB 021211 (Fox 
& Price 2002). The best HETE-2 localization at the time of our discovery announcement 
was 14 arcmin in radius; a subsequent SXC localization produced the 2 arcmin-radius circle 
shown in Figure 1. The position of the optical transient (OT) with respect to the Guide 
Star Catalog II is «=08:08:59.858 and ^=+06:43:37.52 (J2000) with estimated uncertainty 
of 0.1 arcsec in each axis. 

We continued monitoring with P48, the 40-inch telescope at Siding Spring (SSO-40), 
and the Palomar Hale 5-m (P200) telescope with JCAM (Bloom et al. 2003). Photometry 
has been performed using the secondary calibrators of Henden (2002); see Table 1. 

Our optical source was subsequently identified in images obtained by several robotic 
telescopes and other ground-based facilities. In particular, RAPTOR (Wozniak et al. 2002), 
KAIT (Li et al. 2002) and S-LOTIS (Park et al. 2002) detected the source 90 s, 108 s, and 143 s 
after the burst, respectively. In Figure 2 we present the light curve obtained by combining 
data from these and other GCNs with our data. 

Near-infrared - We observed the afterglow in the near-infrared with the Cassegrain 
f/70 low-background secondary and D-78 infrared camera (Soifer et al. 1991) on P200 on 
Dec. 11.55 UT (i^s-band) and with NIRC (Matthews & Soifer 1994) on the Keck-I 10- 
m telescope on Dec. 18.5 UT (JHKs). NIRC observations were photometered relative to 
SJ 9116 and SJ 9134 (Persson et al. 1998), and the P200 zero-point was determined from the 
NIRC observations by reference to two secondary calibrators, SI and S2 (see Fig. 1). Our 
results are given in Table 1. 

Radio - We searched for radio emission with the Very Large Array (VLA)^; our initial 
observation eit t — 2.44 hr represents the earhest radio observation of any GRB to date. We 
used PKS 0734+100 for phase calibration and 3C 147 and 3C 286 for flux cahbration. No 
counterpart to the OT is detected in any single observation. Adding the data between 2002 
December 20 and 2003 January 6, we measure an 8.46-GHz upper limit of 35 /iJy at the 
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position of the OT. Peak flux densities at tfie position of the OT can be found in Table 2. 

Separately, on Dec 12.51 UT we carried out an observation in the 347-GHz band with 
the Submillimetre Common User Bolometer Array (SCUBA) on the James Clark Maxwell 
Telescope (JCMT)^. We used CRL618 as the secondary flux standard. At the position of 
the OT we measure a 347-GHz flux density of 2.0 ± 2.5 mJy. 

3. A Two-Component Light Curve 

To date, only three bursts have had early-time optical detections: GRB 990123 (Akerlof 
et al. 1999), GRB 021004 (Fox et al. 2003) and GRB 021211. GRB 021211 is at any epoch 
one of the faintest detected afterglows; at early epochs it is more than 3 mag fainter than 
GRB 990123, and at later epochs {t > Id) it hes below the great majority of afterglow 
lightcurvcs (Fig. 2). As we show below, this is not the result of extinction in the host galaxy, 
so this burst reminds us that many so-called "dark bursts" may simply be faint (c.f. Fynbo 
et al. 2001; Berger et al. 2002). 

Examining Figure 2, we see that the early behavior of GRB 021211 bears a striking 
resemblance to the light curve of GRB 990123 (c.f. Chornock et al. 2002). The hght curve can 
be considered as composed of two segments: an initial steeply-declining "flash", / oc with 
ttrs ~ 1.6, followed by emission declining as a typical afterglow, ats ~ 1. The corresponding 
terms are "flare" and afterglow for the radio emission (Kulkarni et al. 1999). 

The optical flash is thought to result from a reverse shock propagating into the GRB 
ejecta (Sari & Piran 1999a), while the afterglow arises in the shocked ambient material 
swept up by the ejecta (Sari et al. 1998). For each shock, it is assumed that some fraction 
of the shock energy is partitioned into energetic electrons (eg) and magnetic fields (e^)- The 
electrons are expected to follow a power-law distribution in energy, dN/ d'j gcj'^ for 7 > 7m, 
where 7 is the electron Lorentz factor. Their collective emission follows f„ oc at 
frequencies above the synchrotron peak frequency 1/^ and f^, oc {v/i'mY^^ otherwise. Below 
the self- absorption frequency i^o, synchrotron self- absorption becomes important, and above 
the characteristic cooling frequency i^c, electrons lose their energy rapidly, on the timescale 
of the shock evolution. 

The forward and reverse shocks evolve differently. Once the reverse shock has passed 



^The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the Particle 
Physics and Astronomy Research Council of the United Kingdom, the Netherlands Organization for Scientific 
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through the ejecta shell then no new electrons are accelerated and the shocked electrons cool 
adiabatically: oc t~^^/^^, with the flux at (the "peak" flux) decaying as F,y^ oc ^-^^48 
(Sari & Piran 1999b); for u > the emission decreases even more rapidly. In contrast, 
electrons are continually added to and accelerated at the forward shock; in the case of a 
uniform circumburst medium, F,y^ is constant, and even above Uc the flux declines only slowly. 
Since in a stellar wind-type circumburst medium (p ~ r^^) the reverse shock i/^ is inevitably 
below the optical band, detection of long-lived reverse shock emission for GRB 021211 suffices 
to rule out wind models for this burst (Chevalier & Li 2000). 



3.1. The Reverse Shock Flcire 

The optical flux is already declining at the flrst epoch (Wozniak et al. 2002), when 
t = t* = 90s and Rc = 14 mag, corresponding to / = /* = 7.7 mJy at a mean frequency of 
fo = 4.7 X 10^^ Hz; here the superscript * identifies quantities taken a,t t — 90s. This implies 
that 1/^ < Uo, and constrains the time when the fiare spectral peak i/^ passes through the 
optical bandpass, tm,o: to be tm,o < t*. Since the observed decline a,t t > t* is not very rapid, 
Qirs ^ 2, we have also that u* > u^. 

The constraint on the time of the optical peak flux can be converted into a constraint 
on the relativistic Lorentz factor 7 of the GRB ejecta. Prom Sari & Piran (1999b) we have 
tpeak — 5 JT-o ^^^7300^^-^52^ ^ (thin-shcU case), where tpeak is the time at which the reverse shock 
crosses the GRB ejecta, uq is the circumburst particle density per cubic cm, 7300 is 7/300, 
and -E52 is the isotropic-equivalent burst energy in units of 10^^ erg. Since tpeak < tm,o < t*, 
we find that 7 > 100 Uq^^^ E^i^ . 

Using the prompt emission fluence of 10~^ erg cm~^ (8-40 keV) together with the pro- 
posed redshift of ^ = 1.006 (Vrccswijk et al. 2002), we find an isotropic 7-ray energy of 
Eiso = 3 X 10^^ ergs for GRB 021211. Since this fluence is measured over a relatively narrow 
band, the full /c-corrected fluence (Bloom et al. 2001) may easily reach 10^^ ergs, and it is 
reasonable to assume £^52 ~ 1. 

Following Sari & Piran (1999a), in the absence of synchrotron self-absorption the radio 
flux will follow Fj. oc Jqj, times t < t„jr, and Fj. oc t^"" for times t > t^,!, where tm,r 

is the time when z/„j passes through the radio band, Uj. = 8.5 GHz. The value of tm,r follows 
directly from the equation tj^ j. = tm,o(^^o/i^r)^^^''^, and is equal to 1.35 d in the limiting case 
tm,o = t*. By this time the peak flux will have declined as Fm,T = Fm,oitm,r/tm,o)~^'^^^^ ■ Note 
that we are temporarily ignoring extinction effects, which are discussed in §3.2 below. 

For tm,o < t*, we extrapolate the observed optical flux decay back in time as F^^^ — 
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f*{tm,o/t*) The light curve of the radio flare is thus a direct function of t 
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with the decay past tm,i steepening to Our 3-a upper limits at t = 0.10 d and t = 0.92 d 

are 110/iJy and 66 /^Jy, respectively. For a^s ~ 1.6 (§3.2), the radio flux ai t = 0.92 d is 
predicted to be ~15 /iJy, independent of tm,o- The radio flux at t = 0.10 d is a strong function 
of tm,o, however, and Eq. 1 gives tm,o > 20 s and F^^o < 135 mJy. 

A lower limit on t„i,o implies an upper limit on the Lorentz factor of the ejecta, 7 < 
180 ng -E52 . Following the treatment of Sari & Piran (1999a), however, we find that 
synchrotron self-absorption may be suppressing the radio flux below the level of Eq. 1, 
which would render this limit invalid. 



A variety of afterglow lightcurves are possible (Kobayashi 2000) but the form of sharp 

rise, plateau, and power-law decay is near-universal. Given the bright flare, we cannot 
constrain the rising portion of the afterglow lightcurve. We flt the flare as a power-law decay 
of index a^-s and the afterglow as a constant flux, fm, which undergoes a power-law decay 
with index afg for t > t^, (Fig. 2). 

We flrst examine the constraints on the reverse shock decay. Published GCN data are 
sufficient to constrain a,.^ = 1.63 ± 0.13, and since ctrs = H + i|(^^) (Sari & Piran 1999a), 
this implies that p = 1.85 ± 0.17. 

The gap between the GCN data and our P48 measurements creates fitting degeneracies 
in the full model, so we fix cirs at three successive values, 1.50, 1.63, and 1.76, and make 
fits to the full dataset. Incorporating the range in parameter uncertainties covered in these 
three fits, we find that = 0.80 ± 0.13, = 65 ± 20 /xJy and h = 1700 ± 280 s. Note that 
because the flare dominates the afterglow emission at early times, our fitted should be 
considered a lower limit and our tf, an upper limit. 

We can measure the power-law slope of the spectrum, f3ohs = —0.98, using our B and 
Ks photometry near t = 0.1 d (correcting first for Galactic extinction of Aij=0.07mag; 
Schlegel et al. 1998). If < at this time, then ctfe = 3/4(p - 1) + 1/4, p = 1.7 ± 0.2, 
and P — —p/2 — —0.85 ± 0.1, implying little extinction in the host, Aji — 0.15 ± 0.18 in 
the observer frame. On the other hand, if Uc > i^o at this time, then ctfs = 3/4(p — 1), 
p — 2.1± 0.2, and /3 = — (p — l)/2 = —0.55 ± 0.1, corresponding to an intrinsic reddening of 



3.2. The Forward Shock Afterglow 
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Ar — 0.48 ±0.18 (observer frame). In either case, extinction within the host cannot explain 
the ~2 mag gap between the i?-band hghtcurve and those of most afterglows. 

Since we do not know the location of Vc relative to the optical band, we adopt an average 
extinction correction of — 0.3. This corrects our limiting peak flux of fm ^ 65 /iJy to an 
unextinguished flux of ^ 85 /iJy. 

Prom our model fits, afs = 0.8 ± 0.13, so that Ri 85{tm,o/tb)^'^'^ /^Jy- We can use this 
relation along with Eqs. 4.3 and 4.5 of Sari & Esin (2001) to estimate the circumburst density 
on length scales of r ~ 2jHbC < 10-*^^ cm. We find that no — 0.1e^^^(^^e~^_^(^^ E^2^^^^ cm~^, 
where eg = 0.1ee,-i. Since we have £^52 ~ 1 (§3.1), this implies either a low-density circum- 
burst medium, or abnormally low values of or eg. 

In the absence of self-absorption, the radio afterglow can be expected to peak at time 
t ^ hi'i^o/'i^xY^^'^^ ^ 25 d. In the standard model the synchrotron peak flux remains 
constant as the peak frequency decreases, so we expect a radio peak flux of roughly 85 /iJy. 
Prior to the peak time, Vj- < and F^ ~ t^'^. At the mean epoch of our summed image, 
t = 17.3 d, we expect F,- ~ 70 //Jy, and our 2-a upper limit is F^ < 35 /iJy, a factor of two 
below the predicted level. 

We can estimate the self-absorption frequency (Eq. 4.1 of Sari & Esin 2001) in terms 
of the peak flux, = 0.7 £^52^^^*^ ^q^^^ (Fm/85 fi^jY'^^ GHz. Self-absorption will be effective 
in the case of a high circumburst density - unlikely on the basis of our uq estimate - or if 
Fjn is sufficiently large. In the latter case, at fixed frequency the effects of increasing self- 
absorption win out over the raw increase in flux. Suppressing the radio flux by a factor of 
two requires Fm > 800 /xJy (tm,o ^ 100 s), however, which seems unlikely. 

The discrepancy between the observed and expected fluxes can be due to an early jet 
break, t^et < tm,r (Rhoads 1999; Sari et al. 1999). The optical light-curve, however, shows 
no sharp steepening out to t > lOd (Pig. 2). Moreover, the measured fluence and proposed 
redshift z — 1.006 for GRB 021211 are consistent with minimal beaming for a "standard 
candle" burst (E'lso ~ 5 x 10^° ergs; Prail et al. 2001), so that an early jet- break seems an 
unlikely scenario to explain the suppression of the radio flux. 

Alternatively, the afterglow radiative corrections (Sari 1997; Cohen et al. 1998) may be 
significant. In the slow-cooling regime {vc > I'm), radiative corrections can be substantial if 
ee is large (eg > 0.1), Uc is close to Um, or p is close or equal to 2. We have not derived an 
estimate for eg, but p f« 2 for our models, so this scenario remains an attractive candidate 
for future modeling. 
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Table 1. Optical Photometry for GRB021211 



Date 
(UT) 


AT 
(days) 


Telescope 


Filter 


Magnitude 


Error 








1 V 




n 094 


Dpp 1 1 4Qfi9 






r>a 
1 1 


18 813 


n 04 


r>cin 1 1 /IQQ7 


n 0987 
U.UZo ( 


P4S 


-Da 
i\, 


1 Q DQ'? 


u.uoy 


L^ec ii.ouu^ 


u .uoo^ 


PAS 




1 Q '?4'? 
iy. o^o 


n (171 

U.U ( i 


Dec 11.5120 


0.0410 


P48 


R-^ 


19.286 


0.068 


Dec 11.5139 


0.0429 


P48 


R'^ 


19.479 


0.077 


Dec 11.5218 


0.0508 


P48 


R" 


19.529 


0.082 


Dec 11.5251 


0.0541 


P48 


R" 


19.598 


0.086 


Dec 11.5270 


0.0560 


P48 


R" 


19.755 


0.102 


Dec 11.5332 


0.0622 


P48 


R" 


19.950 


0.110 


Dec 11.5399 


0.0689 


P48 


R" 


20.091 


0.130 


Dec 11.5533 


0.0823 


P200/IRcam 


Ks 


18.01 


0.150 


Dec 11.6025 


0.1310 


SSO40/WFI 


B 


21.877 


0.138 


Dec 11.7083 


0.2370 


SSO40/WF1 


R 


21.096 


0.128 


Dec 12.518 


1.037 


P200/JCAM 


g' 


23.398 


0.047 


Dec 12.518 


1.037 


P200/JCAM 


r' 


23.416 


0.073 


Dec 18.49 


7.019 


K-I/NIRC 


Ks 


22.12 


0.18 


Dec 18.52 


7.049 


K-I/NIRC 


J 


22.77 


0.11 


Dec 18.56 


7.089 


K-I/NIRC 


H 


22.48 


0.16 



^P48 unfiltered observations have been photometered against the 
i?-band cahbration of Henden (2002). 
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Table 2. Radio Observations of GRB 021211 



Date 


AT 


Telescope 


Frequency 


Flux 


Flux Error 


(UT) 


(days) 




(GHz) 


(/^Jy) 


(/^Jy) 


2002 Dec 11.57 


0.10 


VLA 


8.46 


12 


36 


2002 Dec 12.39 


0.92 


VLA 


8.46 


-36 


22 


2002 Dec 13.40 


1.93 


VLA 


8.46 


9 


45 


2002 Dec 13.40 


1.93 


VLA 


22.5 


48 


62 


2002 Dec 15.32 


3.85 


VLA 


8.46 


45 


23 


2002 Dec 16.29 


4.82 


VLA 


8.46 


-6.5 


19 


2002 Dec 20.32 


8.85 


VLA 


8.46 


60 


28 


2002 Dec 22.26 


10.79 


VLA 


8.46 


0.3 


27 


2002 Dec 26.44 


14.97 


VLA 


8.46 


-7.3 


23 


2002 Dec 28.55 


17.08 


VLA 


8.46 


33 


28 


2003 Jan 04.39 


23.92 


VLA 


8.46 


2.4 


24 


2003 Jan 06.28 


25.81 


VLA 


8.46 


46 


24 


Dec 20-Jan 06 




VLA 


8.46 


15 
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Fig. 1. — Discovery image of the optical afterglow of GRB 021211, as obtained with the Near Earth 
Asteroid Tracking (NEAT) camera on the Oschin 48-inch telescope, Palomar Observatory, 21 min 
after the burst (left). The i2-band (f-emulsion) Palomar Digital Sky Survey (DPOSS) image of the 
GRB localization region is shown for comparison (right). The afterglow is marked by two dashes. 
The circle of radius 2arcminutes marks the subsequent SXC localization. The position of SI, one 
of our secondary standards for NIR photometry (§2), is indicated by a carat on the NEAT image; 
S2 is not visible in either image. Coordinates of SI are a=08:09:00.0, (5=+06:43:52 (J2000) and 
coordinates of S2 are a=08:09:00.3, (^=+06:43:38 (J2000). 
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Fig. 2. — Light curve of GRB 021211 and other GRBs, including a two-component flare + afterglow 
fit for the early optical emission from GRB 021211. The data for GRB 021211 are drawn from the 
GCN literature (circles; see text) and this work (squares; Table 1). We fit the data at t < 0.1 d 
only (see §3.2 for details); the dotted line at t > O.ld is merely illustrative. The last two points 
(t > 5 d) are derived from HST observations (Fruchter et al. 2002) and should be largely free of host 
contamination. The data for other GRBs are drawn from the literature (see Berger et al. 2002). 



